Critical behavior of Ginzburg-Landau model coupled to massless Dirac fermions 
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We point out interesting effects of additional massless Dirac fermions with TVf colors upon the 
critical behavior of the Ginzburg-Landau model. For increasing JVf, the model is driven into the 
type II regime of superconductivity. The critical exponents are given as a function of JVf. 
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The critical fluctuations in the Ginzburg-Landau (GL) 
model of superconductors are an old problem in con- 
densed matter physics [Q . While the underlying complex 
order field theory with |0| 4 -interaction is well understood, 
|^| no satisfactory approximation has been found for a 
long time to deal with the additional gauge field. This 
may seem surprising since the Lagrangian is quadratic 
on the gauge field A. One has therefore expected that 
A can be integrated out in a reasonable approximation 
to obtain an effective action with extra terms in the or- 
der field <f>. [EJ This is exactly possible for constant \<f>\ 
where a mean-field approximation for the effective poten- 
tial receives an extra term ~ — \<p\ d in d dimensions. In 
four dimensions, where the model is relevant to particle 
physics, the extra term is ~ |</>| 4 In \4>\ 2 . || Such an ex- 
tra term, if present in the full effective potential, would 
make the second-order phase transition first-order. A 
similar conclusion is derived from a one-loop renormal- 
ization group (RG) calculation in d = 4 — e dimensions, 
which shows no non-trivial charged fixed point [|, || , even 
up to two loops (|, H ■ If the GL model is generalized in 
a such a way as to contain N/2 complex scalars instead 
of one, then non-trivial charged fixed points are found at 
one loop for N > N c — 365.9. Duality arguments, how- 
ever, point out to the existence of a second-order phase 
transition at N — 2 and in the type II regime |Q p"2| . 

A more significant reduction of the critical value of 
N is achieved by a RG approach in a fixed dimension 
d £ (2,4). As we shall see, a one-loop calculation in 
d = 3 reduces N c less than a third of the above value. 
This leads us to expect that non-trivial charged fixed 
points are more accessible in d = 3 than in 4 — e dimen- 
sions. Indeed, this was recently confirmed by the present 
authors || by finding such a fixed point at N = 2 in 
a new three-dimensional RG calculation below T c . The 
success of this approach relies on the explicit presence of 
two mass scales in the problem, defined by the inverse of 
the correlation length £ and penetration depth A. This is 
in contrast to all previous studies which were done in the 
disordered phase at T > T c which has only one physical 
mass scale. It has been proposed to introduce a sec- 
ond scale by assuming different renormalization points 
for each coupling of the theory ||. Such a procedure 
has in fact led to charged fixed points at N = 2, d = 3 
provided the ratio between the two scales is sufficiently 



large ||. However, this ratio must be fixed by conditions 
external to the formalism, and a duality result for a 
tricritical point was used to do so. Thus, in contrast to 
the recent finding for T < T c , the RG-theoretic explana- 
tion of a charged fixed point for T > T c at N = 2 has 
remained obscure. 

The purpose of this paper is to point out the interesting 
modifications of the critical behavior of the GL model 
with N/2 complex field <f> brought about by the presence 
of extra massless Dirac fermion fields ip with Np replica. 
The introduction of massles Dirac fermions is of great 
interest, since effective microscopic models of strongly 
correlated electrons usualy contain them Jl0|, O . 

The bare Lagrangian is assumed to be 
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where the subscript zero denotes bare quantities and 
= d^Al - <9„A°. The labels for the N F fermion 
and N boson replica are omitted. 

The fermions have the effect of modifying the gauge 
field properties of the GL model by giving it an effective 
non-local gradient energy. Indeed, integrating out the 
fermions generates a leading long-wavelength energy 
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In the infrared, this leading term makes the initial 
Maxwell term in (||) irrelevant. Since (Q) gives the gauge 
field a unit dimension instead of the canonical 1/2, the 
charge becomes effectively dimensionless. Thus, by in- 
tegrating out the gauge field for a uniform order field 
</>o = <f>, we obtain the effective potential: 
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where A is an ultraviolet cutoff. Note the important dif- 
ference with respect to the effective potential of the usual 
GL model, where a term \<fi\ 3 is generated in three dimen- 
sions [jlj, giving rise to an apparent first-order transition. 
The limit N F — >■ is singular in this approximation, 
which ignores the Maxwell term controlling the gauge 
field fluctuations for N F = 0. 

For large Np, V c g reduces to the mean-field effec- 
tive potential of the pure |</>| 4 theory. This should be 
not surprising since decoupling by rescaling the charge 
e o - * e o / \/^Vf and taking the large Np limit leads to 
an extreme type II superconductor coinciding with the 
O(N) model. Thus, Np Dirac fermions allows a novel 
interpolation between the usual GL model and the O(N) 
model which runs through different intermediate physical 
systems than the simple limit e\ — * 0. It is therefore an 
interesting problem to study their effect upon the critical 
behavior as the number Np is varied for fixed AT. This 
is what will be done in this paper using RG techniques. 
At one loop and in d — 4 — e dimensions we find that 
for N — 2, which is the physical number for a supercon- 
ductor, an infrared stable charged fixed point exists for 
Np > N-p c = 3.47. We repeat the study at fixed dimen- 
sions d £ (2,4) where we find that the critical number 
of fermions for d — 3 is almost the same: Af c = 4.44 
such that we can give the scheme-independent estimate 
Np c = 4±0.5. Finally, all independent critical exponents 



will be listed as a function of Np . 

Taking into account the Ward identities due to gauge 
invariance, the Lagrangians (§ and (|) can be written in 
terms of renormalized quantities as 
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We define the dimensionless couplings / = /i^Z^e 2 , and 
g = /i~ e Z 2 l Z~ 1 uo, and fix the renormalization constants 
by minimal subtraction of 1/e™ pole terms. The resulting 
one-loop /3-functions are 
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For Np = we recover the usual one-loop /3-functions of 
the GL model [Q, Q]. Note that f3 g is unaffected by the 
fermions, being just the one-loop /3-function of the GL 
model. 

The fixed points lie at 
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where 

A = -2160 - 360N+N 2 + 576N F + 16NN F + 64N F 2 . 

(12) 

Accessible charged fixed points are obtained only if A > 
0. The case of interest for superconductivity is A^ = 2 
for which Eq. ( |l2j ) gives, under the condition A > 0, a 
charged fixed point if 



N F > N F 
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We see that the number of fermions does not need to 
be large in order to produce charged fixed points. A 



schematic flow diagram is shown in Fig. It exhibits 
precisely the fixed point structure expected for the GL 
model j|, ^, Q. It has an infrared stable fixed point 
at (g+, /*), labeled 'SC' in the figure, which governs the 
superconducting phase transition. The zero charge non- 
trivial fixed point labeled 'XY' governs the superfluid 
4 He transition with XY critical exponents. This fixed 
point is unstable for arbitrarily small charge. There is a 
second charged fixed point labeled 'T' which is infrared 
stable only along the line flowing from the Gaussian fixed 
point to it. This fixed point is called the tricritical fixed 
point and the line of infrared stability is a tricritical line. 
The tricritical fixed point has coordinates (g*_, /*). The 
tricritical line separates the left-hand region where the 
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FIG. 1: Schematic flow diagram for the cases where the Np 
and N values are such that charged fixed points are generated. 



phase transition is first-order from the right-hand region 
where it is second-order. The tricritical fixed point in 
the flow diagram is consistent with the proposed phase 
diagram obtained from duality arguments, where the ex- 
istence of a tricritical point was predicted by one of usfl^] 
in 1982, and recently confirmed by Monte Carlo simula- 
tions II . 

Let us compute the anomalous dimension for the case 
Np = 4 and N = 2. This is given by the fixed-point value 
of the RG function 
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Remakably, the anomalous dimension r]^ is negative as in 
the GL model, where it was for a long time a great puzzle, 
explained only recently as a consequence of momentum 
space instabilities in the order field correlation function 

In order to evaluate a second critical exponent such as 
v we need another RG function 



d . (Z m 

At one-loop order, 7 m is found to be 
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The critical exponent v is obtained from the infrared sta- 
ble fixed point value of the function — 1/(2 + j m ) 
Expanded to order e with N = 2 and, say, Np = 4, we 



TABLE I: Critical exponents and values of the Ginzburg pa- 
rameter k at the tricritical and superconducting fixed point 
for e — 1 and N = 2 for several values of TVf . 
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4 


-0.53 


0.506 


1.24/^/2 


1.77/V2 


5 


-0.43 


0.54 


1.09/^ 


2/V2 


6 


-0.36 


0.56 


I.OI/a/2 


2.17/V5 


10 


-0.21 


0.59 


0.82/^ 


2.68/^ 


15 


-0.15 


0.6 


om/s/2 


3.17/V2 


20 


-0.11 


0.61 


0.61^/2 


3.58/^/2 


100 


-0.02 


0.623 


0.29/V2 


7.47/^ 


1000 


-0.002 


0.625 


0.1/V2 


23.19/^ 



have 7^ = 0.024e. In three dimensions, e = 1, such that 
v k, 0.506. 

Below T c , the gauge field acquires a mass whose inverse 
is the penetration depth A. The ratio between the A and 
£ defines the Ginzburg parameter k. Its square can be 
expressed in term of the coupling constants as k 2 = g/2f. 
At the mean-field level type I and type II superconduc- 
tivity are observed for n < 1/V2 and k > 1/V%, rep- 
sectively. Fluctuations will renormalize this separation 
point to k* 2 = g*_/2f*. The value of k 2 at the supercon- 
ducting fixed point is given by n* 2 = g*^/2f*. For N = 2 
and Np — 4, we have 



1.24 

71"' 



1.77 
7T 
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Both values are above the mean- field GL value l/\/2, 
in contrast to the theoretical |l2| and the Monte Carlo 
numbers |l3| in the GL model. 

In Table I we show the values of critical exponents and 
the Ginzburg parameter for N = 2 and e — 1 for grow- 
ing values of N-p. We see that anomalous dimensions 
approach zero as Np becomes large. The anomalous 
dimension tend to zero with increasing N-p much more 
rapidly than the other quantities in Table I. Note that 
K*_ decreases with Np while k+ increases with Np. In- 
terestingly, the critical exponent v does not change very 
much with Np, attaining quickly a limit value v w 0.625. 

The limit Np — > 00 at fixed N can be done analyt- 
ically in the above equations: limjv F ^oo /* \n=2 —> 0, 
limjv P ^oo5-|w=2 -> 0, and Umjv P -«x> 5+|at=2 -* 8ir 2 /5. 
This implies that for many fermions, n* + — > 00 while 
k*_ — > for N = 2. The limiting critical exponents are 
770 = and v = 0.625. Since k* + —* 00 as Np — > 00 at 
N = 2 fixed, this limit is an extreme type II limit in our 
model. The limiting value v — 0.625 equals the one-loop 
result for the XY- model in the e-expansion . 

Let us compare the result in d = 4— e dimensions to the 
fixed dimension RG approach. Instead computing the (3- 
functions for e small, we can set m 2 — and compute 
the Feynman integrals for any dimension d £ (2,4). The 
/3-functions are in this case given at one loop by 
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f = (4 - d){-f + [8N F A(d) + NB(d)]f 2 }, (19) 



TABLE II: Critical exponents and values of the Ginzburg 
parameter at the tricritical and superconducting fixed points 
for d = 3 and N — 2 for several values of Nf- 



P a = (4-d) -g + C(d) 
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5 


-0.36 


0.53 


1.13/V2 


1.6/V2 


6 


-0.31 


0.55 


1/V2 


1.79/V2 


10 


-0.19 


0.59 


0.79/V2 


2.28/V5 


15 


-0.13 


0.6 


0.65/^2 


2.7/V2 


20 


-0.1 


0.61 


0.6/V2 


3.08/^2 


100 


-0.02 


0.623 


0.27/^2 


6.44/^ 


1000 


-0.002 


0.625 


0.09/^2 


20.06/\/2 



A(d) 



T(2-d/2)T 2 {d/2) 
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B(d) = - 



T(l-d/2)T 2 (d/2) 
(4ir) d / 2 r(d) ' 



C(d) 



r(2-rf/2)r 2 (d/2-i) 



(47r)d/2r(d - 2) 
The RG functions 7^ and j m , are at one loop: 

^ = (l-d)(4-d)C(d)f, 
7m = (N+ 2)(d - 4)C(rf) 5 /2 - 70 . 
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Since we are working at the critical point, the RG func- 
tion 7 m above is obtained from an insertion of the com- 
posite field \4>\ 2 into the two-point function. 

As a cross check we set d = 4 — e and expand to first 
order in e, and verify that the /3-functions ([l9|) and (|2fj| ) 
reduce correctly to the previous (||) and , respectively. 

In the absence of fermions, the critical value of N above 
which charged fixed points exist for d = 3 is N c = 103.4, 
much smaller than the value given in the e-expansion, 
N c = 365.9. On the other hand, when we set N = 2 
the critical number of fermions is larger than in the e- 
expansion, being given by Np c = 4.44. On the basis of 
the e-expansion result, we can give the scheme indepen- 
dent estimate as Np c = 4 ± 0.5. 

In Table II we show the values of critical exponents 
and Ginzburg parameter at d = 3 and N — 2 for several 
values of Np. Qualitatively we observe the same behavior 
as in Table I. 

The large Np limit at fixed N = 2 and d — 3 gives 
/* = and limjv F ^oo 9+\n=2 ~ ► 8/5. Thus, we obtain 
the large Np value of v at N = 2 and d = 3: v = 5/8 ~ 
0.625. As before, we have = at Np = 00. 

The fermionic Lagrangian Cp control the influence 
of the thermal magnetic fluctuations, suppressing them 



with increasing Np. An interesting physical case in d = 3 
dimensions has fermion number Np = 10, where the crit- 
ical exponents and values of n at the tricritical and su- 
perconducting fixed points are close to the predicted 
and Monte Carlo -measured values |l3|, [IT], |l8| for the 
pure GL model. The critical exponent v obtained from 
Monte Carlo simulations |17 has the XY model value 



v ~ 0.67, as predicted ]19J from disorder field theory of 
superconductors [p0| . The anomalous dimension 77^ ob- 
tained from Monte Carlo simulations is jig] ~ —0.18. 
The value of n at the tricritical point predicted from the 
disorder field theory is Q k*_ = 0.79/\/2, confirmed by 
recent the Monte Carlo simulations giving k *_ = 0.76/ \f2. 
p3[ The results in Table II, show amazing agreement of 
the theoretical values of 77^ and k*_ in the present model 
with Np = 10. The critical exponent v = 0.59 deviates 
slightly from the one-loop XY result. Presently we don't 
know if this result is just a lucky accident. Anyway, such 
a coincidence deserves further study. 

For A^p < 4 the model has no second-order phase tran- 
sition but interesting physical effects can be expected due 
to chiral symmtery breaking. At d = 3 the Lagrangian 
(|l|) has a chiral symmetry pl| -0 — > exp(i^3^9)tp, with 
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where / is a 2 x 2 unit matrix. In the absence of 
scalar bosons, this symmetry is spontaneously broken 
for A^f < 32/-7T 2 w 3.24 and a fermion mass is dynam- 
ically generated. [^2| According to Kim and Lee [ pr| , 
coupling to bosons reduce this upper bound by factor 
of two and one has Np < 16/ir 2 w 1.62. This lies be- 
low the critical value A^Fc = 4.44 for the existence of 
charged fixed points. The critical behavior described in 
this paper is thus apparently not affected by the chiral 
symmetry breaking. However, we may wonder if the dy- 
namical mass generation in the chirally symmetry-broken 
phase can generate new fixed points in our system. This 
point is very important concerning recent theories of the 
pseudogap state in the cuprate superconductors [p3[ . 
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